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AN IMPROVED LOWER BOUND FOR THE L2-DISCREPANCY 


AICKE HINRICHS AND GERHARD LARCHER 


Abstract. We give an improved lower bound for the L 2 -discrepancy of finite point sets in 
the unit square. 


1 . Introduction and main result 

Let V = {pi,P2, ■ ■ ■ iVn} be a finite point set in the h-dimensional nnit cnbe [0, l)*^. For a 
rectangular box = [0,xi) x ■ • ■ x [0,3;^) in [0,1)'^ let the discrepancy function be dehned 
as 

N 

Dvix) := ^XsAPi) - N -Xi-'-Xd, 

i=l 

where denotes the characteristic function of the set and x := (xi,. .. ,Xd)- So the 
discrepancy function measures the deviation of the number of points of V in the box B^ 
from the fair number of points N ■ Xi ■ ■ ■ Xd ='■ N\Bx\. 

The L2-discrepancy of V is the L2-norm of D-p{x), given by 

\\Dr\L 2 \\ = ([ D^x)dx 

By a celebrated result of K.F. Roth (see also |KNj ) it is known, that for every dimension 
d there are constants 0 < < c'd such that 

(i) for all A-point sets V in [0,1)'^ we have 

||Up|L2il >Q-(logA)^ 

and 

(ii) for any N >2, there exists an A-point set V in [0,1)'^ such that 

||L)p|L2|| < c'd • (log A)~ . 

Hence, it makes sense to ask for the positive real q which is dehned by 

Cd'.= inf inf IID75IL2II ■ (log A) 2 

N>2#r=N 

where the second inhmum is taken over all A-point sets V in [0,1)'^. 

The second author is partially supported by the Austrian Science Fund (EWE), Project F5507-N26, which 
is a part of the Special Research Program “Quasi-Monte Carlo Methods: Theory and Applications”. 

1 






2 


AICKE HINRICHS AND GERHARD LARCHER 


In this work we are mainly interested in dimension 2, i.e., in C 2 - To discuss known results, 
let us also introduce the constants 

_ ^-1 

ttd ■= liminf inf ||Zi)-p|L 2 || ■ (logN) ^ 

N^oo #P=N 

and 

d-1 

bd := limsup inf ||Zi) 75 |L 2 || ■ (logiV) ^ ^ 

AT^oo #'P=N 

which are the interesting constants if we want to study the asymptotic behavior of the 
L 2 -discrepancy for iV —)■ oo. Obviously, 

C-d ^ ^ bd- 

The best constructions known so far yield the upper bound 

02 < 0.17907... 

This is from a recent construction in |FPPSj using generalized scrambled Hammersley point 
sets. There is numerical evidence obtained in jBTRj that the Fibonacci lattices are slightly 
better and give 

02 < 0.176006... 

But until now this is not proved. 

The best lower bounds known so far are 


( 1 ) 

0.038925 ... < 62 

and 


(2) 

0.03276 ... < C2. 


These bounds are shown in |HM] . Actually, there the hrst lower bound was even claimed 
to hold for C 2 . However, there is a small inaccuracy in the proof of the lower bound in jHM] : 
At the end of Section 3, when the optimal value for the parameter 7 is determined, in the 
paper the maximum of the function — 2~^7~^y^ for 1/2 < y < 1 was determined. In 

fact, however, one does not need the maximum but the minimum of this expression, which is 
attained for y = 1. If we carry out the calculation in the correct way, we obtain the correct 
estimate (|2|). But taking the maximum still gives the lower bound (p!]). The estimates in 
[HM] for d > 2 have also to be adjusted in the same way. 

These lower bounds are also valid for the weighted L 2 -discrepancy. In this paper we 
concentrate on the case of equal weights. Here our aim is to improve these lower estimates. 
We show 

Theorem 1. 

C 2 > 0.0515599 ... 

That means: For every N-point set V in [0,1)^ we have 

\\Dr\L2\\ > 0.0515599 ... ■ ^log N. 













AN IMPROVED LOWER BOUND FOR THE L 2 -DISCREPANCY 


3 


Moreover, 

h > 0.0610739.... 

The proof is a variant of the proof given in |HM] but we go one step deeper in the analysis. 
The paper is organized as follows. In Section we collect the necessary facts on Haar- 
coefficients of the discrepancy function. In Section we describe the steps of the proof 
omitting some technical details. These are given in Section 


2. Haar-coefficients of the discrepancy function 


For j = (ji,j 2 ) e { — 1,0,1,2,.. - Y let be the set of m = (mi, m 2 ) G (0,1,.. - Y with 
0 < m, < 2 ^ and let |j| = max(0, ji) + max(0, j 2 )- For such j and m, we consider dyadic 
intervals 

m 2 m2 + l\ 

2^2 ’ 2 I 2 / ’ 

We denote by 1^® left-lower, by 1^® left-upper, by I~Y the right-lower, and by 
I~Y the right-upper quarter of Ij^m- These are again dyadic intervals with a quarter of the 
area of Ij^m- The Haar-functions hj^m are supported by Ij^m and satisfy = 1 on 
and 1 hj,m = — 1 on and ■ Note that all dyadic intervals with hxed j are 
congruent. Therefore we call j the shape of the interval Ij,m- 
Then by Parseval’s equation we have 

( 3 ) = 5^231 5 ^ 

j mGBj 


mi mi -|- 

2 F ’ 2 F ) ^ 


where 


f^j,m / D'p(yX) ■ hj ,m(^Y 
[ 0 , 1)2 

are the Haar-coefficients of the discrepancy function. 

In |HMj the above sum was estimated from below by summation over all j and m such 
that Ij^rn contains no point of V. For these j and m the Haar-coefficient fij^m has a very 
simple form. In this paper we also consider j and m such that contains exactly one 
point of V, and so we will obtain an improved lower bound. 

The following two lemmas were already stated in [Hj : 


Lemma 2. For j G Nq and m G Dj we have 



X1X2 ■ hj^rn{,Xi, X2) dxi dX2 


2-2|i|-4 


Lemma 3. For z = (zi, Z 2 ) G [0, 1)^, let Cz := [zi, 1) x [ 2 : 2 ,1) and let xc^ be the characteristic 
function of Cz- Let j G Ng and m G he such that z ^ Ij^m- Then 



XcAY ■ hj,m(T) da; 


0 . 










4 


AICKE HINRICHS AND GERHARD LARCHER 


In the following we will also have to consider the case when ^ G Ij^m- 
Lemma 4. With the notation of Lemma^ let now z = (zi, Z 2 ) E Ij,m- Then 

- if) (J2 - i») ilzeO;:, 
h - if) (=ig^ - ^2) i!z(=W 


Tj,m ■ 


'[ 0 , 1)2 


XcAx) ■ hj^rnix) dx = 


j,m 5 


(=a+i - 21) (22 - if) ifzeq;:, 
- 2,) (!!|±i - 22) 2/2 6 /^-. 


The lemma is proved by simple calculations, see also Lemma 3.3 in 0. As a corollary of 
the three lemmas we obtain immediately: 


Lemma 5. For the Haar-coefficients of the discrepancy function Dp we have 

t^j,m 


—N2 4 if contains no point from V, 

h'pm ~ if Ij,m contains exactly one point z from V. 


3. Main steps of the proof 

In this section, we describe the high level structure of the proof. The details will then be 
presented in the next section. 

Let N be the number of points in V. Let M G No and 0 < k < 1 be such that N = 2^+'". 
For i G Ng, let |j| = ji + 72 be the level of the dyadic intervals For r G No, let 

^r(j) = {(j, m) : mE and (IV = r} and ar(j) = ifA{j) 

be the set of indices of intervals of shape j containing exactly r points of V and its cardinality, 
respectively. Observe that we have 21-^1 intervals for hxed j and altogether N = 2^+'^ points, 
so 

00 00 

(4) and ^ ra^(j) = 2^+^. 

r=0 r=0 

Since the intervals of a given shape j are mutually disjoint, this implies 

(5) ao(j) > 2l^'l -N = 2l^'l - 2^+^. 

Finally, for i E Nq, let 

= IJ Ar{j) and afT) = ifAfT) = afj)- 
lil=^ \3\=^ 

Since there are £ + 1 different shapes j with level |jj = £ + 1, we obtain from Q 

00 CO 

= {i + 1)2^ and rar{i) = {i + 1)2^"*"^. 

r=0 


( 6 ) 
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We begin with recalling the main idea of the proof of the lower bound in |HM] . Using in 
Parseval’s equation ([^ only the intervals containing no point of V together with Lemma 
and (|^ we have 

CO CO 

(7) \\Dv\L2\\^ >5^2^ (^ + 1)2^ (2^ - 2^+^) 


This sum can then be evaluated and estimated. 

To improve upon this estimate, it is not enough to consider further intervals separately 
since each Haar coefficient of an interval containing points of V can be zero. So we have 
to bundle together some Haar coefficients. In our approach, we consider intervals Ij^m of 
level |j| = M or |j| = M + 1 containing one point together with the two intervals of level 
M + 1 and M + 2, respectively, that are contained in and contain the same point. More 
formally, for {j,m) G Hi(£), dehne 


( 8 ) 




_ 2 I 2 

f^j,m ' 




+ /i 


j" ,m" 


where f and j” with |j'| = |j"| =1+1 are distinct, C Jym and Ij^m HV = 

r\V = n V. For u = 0,1,2, let us call an interval of level £ + 1 with 

G Hi(£ + 1 ) of type u, if exactly u of the two intervals Ij^m of level i containing it 
satisfy {j, m) G Ai{j). That means that exactly 2 — m of these two intervals contain at least 
one additional point of V. We also define for ^>l 


G Ai{i) : Ij^rn is of type u} and 6 „(£) = #S„(£). 
It follows from these definitions that 


cii(£ + 1 ) — + 1 ) + bi[i + 1 ) + b 2 {i + 1 ) and 2 ai(£) — bi(^i + 1 ) + 2 & 2 (f + !)• 

for i G Nq. Indeed, the hrst identity counts the intervals in level £+1 containing exactly one 
point of V in two different ways. For the second identity, consider the bipartite graph with 
intervals in level £ and £ +1, respectively, containing exactly one point of V as the two sets 
of vertices and draw an edge between an interval of level £ and an interval of level 
£+1 il Ij\m' is contained in Ij m- Then the degree of each interval of level £ is exactly 2 
and the degree of each interval //,m' of level £ + 1 is m if (j', m') is of type u. Counting the 
edges in two different ways yields the second identity. This further implies 

(9) 26o(£ + 1) + 61 (£ + 1) = 2ai(f + 1) — 2ai(£). 

The crucial improvement of the estimate ([^ is now 

00 2 

(10) ||C„|L2f>^2' dm+ E + ^ (2 -m)2"pp„. 

e =0 ij,m)£Ao(l) (i,«i)sAi(M) 11=0 (j,m)eS„(M+l) 

Considering ([^, we see that this inequality indeed holds since each Haar coefficient of an 
interval in level M gets at most a weight 2^, each Haar coefficient of an interval in level 
M + 2 gets at most a weight and, by dehnition of Bu{M + 1), each Haar coefficient in 
level M + 1 gets at most a weight 2 ■ 2^ = 2 ^+^. 
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To estimate the pj^m we find in the next section an explicit function 7 : [—1,1] —)■ ]R_|_ such 
that the following lemma holds. 

Lemma 6 . We have > 7(«) if {j,rn) G Ai{M) and 2‘^^pj^rn > 2~‘^'^{k — 1) if 

(j,m) e Ai(M + 1 ). 


\Dv\L 2 f >2 -^T.i + T.2 


Now it follows from (10), ([^, Lemma Lemma and ([^ that 

( 11 ) 

with 

( 12 ) 

and 
(13) 


El = ao{M)2^^-^ + ao{M + 1 ) 2 ^'^-^^ + ai{M)j{K) 

+ + 1 ) + 26o(Tf + 1))2 ^ 7 (^ — 1 ) 

= ao(M)22"-8 + ai(M)(7(/s)) - 2-^-f{K - 1)) 

+ ao(M + 1 ) 2 ^'^-^^ + ai(M + 1 ) ■ 2-^'y{K - 1 ) 


S 2 = ^ 2 ^(£+ 1 )( 2 ^- 2^+'')2 


2M+2k-4^-8 


e=M +2 


Now the next Lemma shows that, under the conditions in ([^, Ei is minimized if the points 
are as evenly distributed in the boxes of level M and M + 1 as possible. 


Lemma 7. For r G No and i = M, M + 1, let ar{i) G No be such that the conditions in (§ 
are satisfied. Then we have 

ao(M)22-« + ai(M)(7(K)-2-i7(ft;-l)) > (M + 1)2'^0 m(/^) 

ao(M + l)22"-ii + ai(M + l)-2-S(«-l) > (M + 2)2'^0 m+i(«:) 

with 


0 m(«:) = (2 - 2 ") ■ (7(«) - 2 -S(« - 1 )) 

= ( 2 - 2 ")- 22 "-ii + 2 "- 2 -S(«- 1 ). 


Now (10) and Lemma [pimply 


(14) 


Dr\L2f > S; + S '2 


with 

(15) s; = (M + 1)(2 - 2") ■ (7(/s:) - 2-S(« - 1)) + (M + 2)2" ■ 2-S(« - 1) 
and 

00 

S 2 = 2^(^ + 1)(2^ “ ]\ i ' j 2 ‘^ m + 2 k - u -8 > (jvf + 2) (3"^ ■ 2^""^ - 7 -^ ■ 2^""®) . 

e=M+i 


(16) 
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The last inequality is the direct computation done in |HMj . Observe that the summation 
in S 2 starts at £ = M + 1 instead of £ = M + 2 for S 2 . We now integrated the summand 
corresponding to the intervals not containing points of V into this sum. This has the ad¬ 
vantage that this sum is exactly the same sum already estimated in |HMj . The complete 
improvement is now contained in 

Altogether, we obtain the improved bound 

111^1^117211 > (in -|- l)A(/i;) > log iV ■ -—^ 

log 2 

with 


(17) A(k) = 3 -^ ■ 22^-® - 7 -^ ■ 2^'^-® + (2 - + (2^^ - 1)7(k - 1). 

A hnal analysis shows that 


(18) 


inf 

Ke[o,i) 



0.0515599... 


and 


sup 

Ke[o,i) 



and hnishes the proof of Theorem [T] 


0.0610739... 


4. Details of the proof 

In this section, we prove Lemmas and and provide the analysis showing (18). 

Proof of Lemma\^ Let (j, m) G Ai{M) and assume that the point p = {pi,P 2 ) of V which is 
in is contained in the left lower part of Ij^m, i-e., with j' = (ji -|- 1 , j 2 ) and j" = {ji,j 2 + 1 ) 
we have m' = ( 2 mi,m 2 ) and m" = (mi, 2 m 2 ). The three other cases are treated in exactly 
the same way. 

This means that 

2^m+k 

Pj,m = xy- 227 + 2 i 2+4 ’ 

where x := pi — ^ and y ■= P 2 — ^ and hence 0 < a; < and 0 < y < ^ 3 ^. 

We have now to distinguish between four cases: 

1 


Case 1: 
Case 2: 
Case 3: 
Case 4: 


271+2 

1 


0 < a: < 
0 < a: < 

< X < 


< X < 


271+2 ’ 
1 

271+2 ’ 
1 

271+1 ’ 
1 


272+2 


0<y < 


<y < 


271+2 — “ ' 2ii+i ’ 

Case 1: In this case, by Lemma we have 

/ 2 ^ m+k \ 2 / 2M-1-K \ ■ 

- i^y ~ 2231+232+4 j +‘^i^y- 22ii+2j2+6 j 


= 2 : 


272+2 


2M+4 —K 


0<y < 


<y < 


272 +2 
1 

272+1 

1 

272 +2 
1 

272+1 


+ 2{z 


2M+6-K 
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with Q < z ■.= xy < Simple analysis shows that pj^m attains its minimnm ^ 2 mAi- 2 k 

for z = 

Case 2: In this case, again by Lemma we have 

~ ~ ~ ^ 22 ^ 2+1 ~ y^ ~ 2 ^+ 6 -«:^ ' 

Hence 

(ph+j2+2\^ _ c^2M+4: 

= +(«( 1 -/^)-^) =:/(«>/^) 

with 0 < a := 2P+^x < f and f < /3 := 2 ^ 2 + 1 ?/ < 1. 

Simple analysis of /(a, /3) shows that / in this region for a and [3 has minimal valne ^^p 
for a = ^ ■ 2^" and /9 = |. Hence, in Case 2, we always have 

1 9-22^ 

Pj,m ^ 22M+4 g]^2 

Case 3: This case is treated in the same way as Case 2 and also resnlts in 

1 9 ■ 2^'' 

Pxm > 

Case 4: In this case, again by Lemma we have 

Pp'^ ~ (^^y ~ 2M+4-k^ ^2^2+1 ~ y^ ~ 2M+6-k^ (^(^2T+T ~ y ~ 2M+6 -k^ ■ 

Hence 

= (a/? - 2^^ + (^a(l -/3) - + (^(1 - a)/3 - =:f7(a,/5), 

now with 

- < a := 2P~^^x < 1 and - < B := 2^‘^^^y < 1. 

The analysis of g{a,(3), i.e., the determination of mma,i 3 gic(, (3) is explicitly possible, bnt 
leads to an eqnation of degree 4. The analysis was carried ont with the help of Mathematica. 
We will not give the details of this in principle elementary bnt elaborate analysis. 

It turns out that g attains its minimum for each k, with 0 < a < 1 for some a = (3 = /3 {k). 
The function /3{k) can be explicitly expressed with fourth-roots. Hence, 

h{K) := mmg{a,/3) = g{/3{K),/3 {k)) 

aP 

for each k is an explicitly given (but ugly) function. Its form is given by the dashed line in 
Figure 1. 

























AN IMPROVED LOWER BOUND FOR THE L 2 -DISCREPANCY 


Hence, in Case 4, 

Pj,m > 22M+4 ■ ^(^)- 
So, altogether we have in any case 

Pj,m > min (3 ■ 9 ■ 2~'^h{K)) = min (9 ■ 2“^h(K)) =: 7 (k). 

So, 7 (k) is the minimnm of the two functions 2~‘^h{K) (dashed line) and 9 ■ (undashed 

line) in Figure 1 and it is explicitly given. 



Figure 1. The functions 2 ^/i(k) and 9 ■ 2“^^ 

The inequality for {j,m) G Ai{M + 1) immediately follows by noting that we can proceed 
in the same way as before by replacing nhj k — 1 and M by M + 1 now. □ 

Proof of Lemma\^ For r G No, let > 0 be such that 

(19) ^^ar = l and ''^^ror = a. 

r>0 r>0 

Furthermore, let a,/3 > 0 be such that a > 2f3. Then we have 

(20) aao + /3ai > a{l — a) + /3a if 0 < a < 1 

(21) aoo + /3ai > /3{2 — a) if 1 < a < 2. 


Indeed, (20) follows directly from the identity 

aoo + f3ai = a{l — a) + 13a + (a(r — 1) — (3r) 


r>2 
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and a{r — 1) — /3r > /3{r — 2) > 0 for r >2. Similarly, (20) follows from the identity 

2ao + oi = 2 — (T + — 2)ar 

r>3 

together with aa^ + (3ai > (3{2aQ + Oi). 

Now the ineqnalities to be proved follow from (20) and (21) by considering the special 
cases 


ar = {{M + 1)2^) ar{M), a = 2^^ e [1, 2], a = 2^^-®, (3 = 7(k) - 2-S(k - 1) 


and 


ar = {{M + 2)2^+^) \r{M + 1), a = 2^“^ e [0,1], a = ^ = 2-^j{k - 1). 

Indeed, the eqnalities (IT^ follow from Moreover, in the first case 


2(3 = 2-f{K) - 7 (a - 1) < 27(k) < 9 ■ 2^^-^^ < 2^^-® = a 


and in the second case 


2/1 = 7 (k - 1) < 9 ■ 2'^^-^^ < 2^''-^^ = a 


also hold. Finally, /I > 0 is obvious in the second case and is easily checked in the hrst 
case. □ 



Figure 2. The function A(k) 


Proof of (18). It is an elementary task to analyze the function A(fi:) on the unit interval. 
With the help of Mathematica and basic differential calculus we see that there is a Aq ^ (0) 1) 
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with A monotonically increasing in [0, kq) and monotonically decreasing in (kq, !]• Moreover, 
we have 

SI 7 

A(0) = A(l) = 3-^ ■ 2-^ - 7-^ ■ 2-^ + 7(0) = 

Hence, 

317 

niin A(«) = min (A(0), A(l)) = —— = 0.00184268 ... 

Ke[o,i] 112032 

(see the illustration of A(k) in Figure 2) and consequently 


/A(/t;) 

inf = 0.0515599... 

Ke[o,i) y log 2 


Moreover, numerical computation gives kq 
hnally leads to 


sup 

Ke[o,i) 



= 0.5705243... and A(ko) 



0.0610739... 


0.00258545 ... which 


□ 


We close with some remarks and open problems. The lower bounds in |HMj are valid for 
the weighted discrepancy. This is mainly due to the fact that only dyadic boxes are used 
that do not contain points of the pointset. To show the better bounds in this paper we need 
to consider also dyadic boxes containing points. This does not seem to easily generalize to 
the case of the weighted discrepancy. Also, to extend the approach to dimension d > 2 seems 
to be technically difficult since there are more cases to consider than in the proof for d = 2. 
It is also an open problem if possibly dd = bd, i-e. if the limit 

lim inf IlD-plLoll ■ (logA)~^~ 

N^oc#r=N I II V / 

exists. 
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